TRI-PP-95-86 
E^Tth/9511021 

Recovering Relativistic Nuclear Phenomenology from the 
Quark-Meson Coupling Model 

Xuemin Jin and B.K. Jennings 
TRIUMF, 4004 Wesbrook Mall, 
Vancouver, British Columbia, Canada V6T 2 A3 

Abstract 

The quark-meson coupling (QMC) model for nuclear matter, which describes 
nuclear matter as non-overlapping MIT bags bound by the self-consistent 
exchange of scalar and vector mesons is modified by the introduction of a 
density dependent bag constant. The density dependence of the bag constant 
is related to that of the in-medium effective nucleon mass through a scaling 
ansatz suggested by partial chiral symmetry restoration in nuclear matter. 
This modification overcomes drawbacks of the QMC model and leads to the 
recovery of the essential features of relativistic nuclear phenomenology. This 
suggests that the modification of the bag constant in the nuclear medium may 
play an important role in low- and medium-energy nuclear physics. 
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Since quantum chromodynamics (QCD) is believed to be the correct theory underlying 
strong interactions, the physics of nuclei is, in essence, an exercise in applied QCD. Building 
connections between observed nuclear phenomena and the interactions and symmetries of 
the underlying quark and gluon degrees of freedom has become one of the principal goals 
of nuclear theorists. However, due to the complexity of the low-energy QCD, knowledge of 
QCD has had very little impact, to date, on the study of low- and medium-energy nuclear 
phenomena. 

A reasonable consensus is that the relevant degrees of freedom for nuclear physics at low 
energy scales are hadrons instead of quarks and gluons. One general approach, relativistic 
nuclear phenomenology which has gained much credibility during last twenty years, is to 
treat nucleons in nuclear environment as point-like Dirac particles interacting with large 
canceling scalar and vector potentials. This approach has been successful in describing the 
spin-observables of nucleon-nucleus scattering in the context of relativistic optical potentials 
0||. Moreover, such potentials can be derived from the relativistic impulse approximation 
0. The relativistic field-theoretical models based on nucleons and mesons, QHD, also feature 
Dirac nucleons interacting through the exchange of scalar and vector mesons ||. QHD, at 
the mean-field level, has proven to be a powerful tool for describing the bulk properties of 
nuclear matter and spin-orbit splittings of finite nuclei 0. It is known that the large and 
canceling scalar and vector potentials are central to the success of the relativistic nuclear 
phenomenology. Recent progress in understanding the origin of these large potentials for 
propagating nucleons in nuclear matter has been made via the analysis of the finite-density 
QCD sum rules ||. 

Given the wide success of the Dirac approach in describing various low- and medium- 
energy nuclear phenomena and the support from the finite-density QCD sum rules, it is a 
challenge to study the relevance of the quark structure of the nucleon to the dynamics of 
normal nuclei. A few years ago, Guichon [H proposed a quark-meson coupling (QMC) model 
to investigate the direct "quark effects" in nuclei. In this model, nuclear matter consists of 
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non-overlapping MIT bags interacting through the self-consistent exchange of mesons in the 
mean-field approximation, and the mesons are directly coupled to the quarks. This simple 
QMC model has been refined by including nucleon Fermi motion and the center-of-mass 
corrections to the bag energy || applied to variety of problems |7Hl"3"11 . (There have been 
several works that also discuss the quark effects in nuclei, based on other effective models 
for the nucleon HT4|). 

Although it provides a simple and attractive framework to incorporate the quark struc- 
ture of the nucleon in the study of nuclear phenomena, the QMC model has a serious short- 
coming. It predicts much smaller scalar and vector potentials for the nucleon than obtained 
in relativistic nuclear phenomenology and finite-density QCD sum rules. Unless there is a 
large isoscalar anomalous coupling (ruled out by other considerations) this implies a much 
smaller nucleon spin-orbit force in finite nuclei. To lowest order in the nucleon velocity and 
potential depth the nucleon spin-orbit potential can be obtained in a model independent 
way from the strengths of the scalar and vector potentials. The spin-orbit potential from 
the QMC model is too weak to successfully explain spin-orbit splittings in finite nuclei and 
the spin-observables in nucleon-nucleus scattering. 

We observe that the bag constant is held to be at its free space value in the QMC model 
for nuclear matter. This assumption can be questioned. In the MIT bag model, the bag 
constant denotes the vacuum energy (relative to the perturbative vacuum), which contributes 
~ 300 MeV to the nucleon energy and provides the necessary pressure to confine the quarks. 
Thus, the bag constant is an inseparable ingredient of the bag picture of a nucleon. When a 
nucleon bag is put into the nuclear medium, the bag as a whole reacts to the environment. 
As a result, the bag constant may be modified. 

There is little doubt that at sufficiently high densities, the bag constant is eventually 
melted away and quarks and gluons become the appropriate degrees of freedom. It, thus, 
seems reasonable that the bag constant be modified and decrease as density increases. More- 



over, it is argued in Ref. [15] that the MIT bag constant is related to the energy associated 
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with the chiral symmetry restoration (the vacuum energy difference between the chiral- 
symmetry- restored vacuum inside the bag and the broken phase outside). Since chiral sym- 
metry is partially restored in nuclear medium |L£j , the in- medium bag constant should drop 



relative to its free space value 15|,16|1. This physics has been ignored in the QMC model. 



In this Letter, we shall introduce the in-medium modification of the bag constant in 
the QMC model for nuclear matter. We relate the in-medium bag constant to the in- 
medium effective nucleon mass through a scaling ansatz suggested by partial chiral symmetry 
restoration in nuclear matter. We find that the essential features of the relativistic nuclear 
phenomenology, in particular the large canceling scalar and vector potentials and hence 
strong spin-orbit force for the nucleon, can be recovered when the decrease of the bag 
constant in medium is taken into account. This suggests that the drop of the bag constant 
in nuclear medium relative to its free space value may play an important role in low- and 
medium-energy nuclear physics. 

In the QMC model, the nucleon in nuclear medium is assumed to be a static spherical 
MIT bag in which quarks interact with the scalar and vector fields, a and ZJ, and these 
fields are treated as classical fields in the mean field approximation. (Here we only consider 
up and down quarks.) The quark field, ip q (x), inside the bag then satisfies the equation of 
motion: 



iP q (x) = , (1) 



where m° q is the current quark mass, and g q a and g% denote the quark-meson coupling con- 
stants. The normalized ground state for a quark in the bag is given by 



jo(xr/R) ^ 
K i(3 q a- rjx{xr/R) ) 



Xq (2) 



47T ' 



where e q = Vt q + g^u R and (3 q = y(Q q — Rm*)/(Q q + Rm*), with Vt q = ^ x 2 + (Rm*) 2 , 
m* = m' — g^a, R the bag radius, and \ q the quark spinor. The normalization factor is 
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given by Af~ 2 = 2R 3 Jq{x) VL q {VL q — 1) + Rm*/2 /x 2 . The x value is determined by the 
boundary condition at the bag surface, jo(x) = j3 q ji(x). 

The energy of a static bag consisting of three ground state quarks can be expressed as 

E hag = 3^-^ + ~7rR 3 B , (3) 

where Z is a parameter which accounts for zero-point motion and B is the bag constant. In 
the calculations to follow, we use Ro, B and Z to denote the corresponding bag parameters 
for the free nucleon. After the corrections of spurious cm. motion in the bag, the effective 
mass of a nucleon bag at rest is taken to be 



M* N = y/^ag - (PL) , (4) 

where (p 2 m ) = X^(Pg) an d (p 2 q ) is the expectation value of the quark momentum squared, 
(x/R) 2 . 

The QMC model assumes that both Z and B are independent of density. (The bag 
radius is determined by the equilibrium condition for the bag, see below.) This assumption is 
unjustified. As argued in Ref. ||15|| , the MIT bag constant is essentially the energy associated 



with the chiral symmetry restoration. Therefore, one expects the bag constant to drop 
relative to its free space value as a consequence of partial chiral symmetry restoration in 
the nuclear medium. Of course, one has to invoke model descriptions in order to obtain a 
quantitative estimate for the reduction of the bag constant. According to the scaling ansatz 
advocated by Brown and Rho W7\, the bag constant should scale like B/Bq ~ $ 4 |T3|,|TB 



Here $ denotes the universal scaling, and $ ~ m*/m p ~ /*/ '/„. ~ [M^ / M^) 2 ^ has been 
suggested in Ref. fT5|]. Motivated by these suggestions, we introduce the following scaling 



ansatz 

B /M* N \ K/3 



(5) 



B \M N 

for the in-medium bag constant. The case k = corresponds to the usual QMC model. 
Note that combining Eqs. (J§), and (El) yields a self-consistency condition for B. In 



principle, the parameter Z may also be modified in the nuclear medium. However, unlike 
the bag constant, it is unclear how Z changes with the density as Z is not directly related 
to chiral symmetry. Here we assume that the medium modification of Z is small at low and 
moderate densities and take Z = Zq. 

The bag radius is determined by the equilibrium condition for the bag, dM^/dR = 0. In 
free space, one may fix at its experimental value 939 MeV and use the equilibrium condi- 
tion to determine the bag parameters. For several choice of bag radius, Rq = 0.6, 0.8, 1.0 fm, 
the results for S 1/4 and Z are 188.1, 157.5, 136.3 MeV and 2.030, 1.628, 1.153, respectively. 

The total energy per nucleon at finite density, p^, can be written as M 
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tot 



J** #kyfM$ + W + £p Pn + 



m 



°a 2 



(6) 



(2tt) 3 P W v " ' 2ml r " ' 2 PN 

where 7 is the spin-isospin degeneracy, and 7 = 4 for symmetric nuclear matter and 7 = 2 for 
neutron matter. Here we have used that the mean field u; created by uniformly distributed 
nucleons is determined by baryon number conservation to be 0-0] 

3 ql Pn 9u Pn 



mi 



(7) 



where g u = 3g*. The scalar mean field is determined by the thermodynamic condition 
(d E tot /da) RpN = 0, which yields the self-consistency condition 



g 2 ^ 

m* [I-k) 6 



d 6 k 



M* N 



M* N 2 + k 2 



(8) 



where g a = 3 g% and 

-^bag 



C(a) 



M* N 



1 - 



-^bag R . 



S(a) + 



m' 



£1 



bag 



1 - 



K E h 



3 Mff 3 



ttR 3 B 



-1 



(9) 



with 



n g /2 + Rm*(Q q -l) 
n q (n q - 1) + Rm* /2 



(10) 



We now turn to present numerical results. For simplicity, we take m° q = 0. The coupling 
constants g a and g u are chosen to fit the nuclear matter binding energy (—16 MeV) at 



TABLE I. Coupling constants and nuclear incompressibility K (in MeV) at various n values. 
The case of k = corresponds to the simple QMC model and the last raw gives the result of 
QHD-I. Here we have used m a = 550 MeV and = 783 MeV. 



k value 




i? =0.6 fm 






iio=0.8 fm 






#0=1.0 fm 






gll^ 


gll^ 


K 




gll^ 


K 


9S/47T 


gll^ 


K 


K = 


20.17 


1.56 


223 


21.87 


1.14 


200 


22.47 


0.96 


189 


K = 7 


4.40 


4.41 


356 


4.51 


3.51 


314 


4.54 


3.18 


304 


K = 8 


3.13 


5.72 


426 


3.14 


4.50 


359 


3.14 


4.07 


350 


K = 9 


2.10 


8.32 


650 


2.04 


6.33 


471 


2.00 


5.67 


415 


K = 10 


1.21 


14.89 


2233 


1.20 


11.17 


1058 


1.16 


9.74 


918 


QHD-I 


8.45 


12.84 


540 


8.45 


12.84 


540 


8.45 


12.84 


540 



the saturation density (p% =0.17 fm™ 3 ). The resulting coupling constants and the nuclear 
incompressibility are listed in Table 1. We note that while the scalar coupling decreases, the 
vector coupling and the nuclear incompressibility increase as k increases. 

This result can be understood from the scaling ansatz Eq. (|5|). When k = 0, the strength 
of the vector field is much smaller than that required in QHD-I. This, in Refs. ||[7|], is at- 
tributed to the repulsion provided by the cm. corrections to the bag energy. When k > 0, 
Eq. (|5D provides a new source of attraction as it effectively reduces M^. Consequently, 
additional vector field strength is required to balance the nuclear matter. The decrease of 
the scalar coupling with increasing k is due to the increasingly strong attraction from the 
dropping bag constant. (When k > 10, the self-consistent solution around pjq = p° N does 
not exist as the attraction gets too strong.) The rapid increase of the nuclear incompress- 
ibility with increasing k results from that the contribution of the vector field to the nuclear 
incompressibility is proportional to g^. 

The total energy per nucleon for symmetric nuclear matter is presented in Fig. [1] for 
various k values, with Rq = 0.8 fm. The result from QHD-I is also plotted for comparison. 
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FIG. 1. Energy per nucleoli for symmetric nuclear matter as a function of the medium density, 
with Rq = 0.8 fm. The five curves correspond to k = (solid), 7 (long-dashed), 8 (dott-dashed), 9 
(short-dashed), and 10 (dotted), respectively. The result from QHD-I is given by the open circles. 

The usual QMC model (k = 0) predicts a much softer equation of state for the nuclear 
matter than in QHD-I. As k gets larger, the equation of state becomes stiffer, and when 
k ~ 9.5 the equation of state is essentially the same as the one predicted in QHD-I. The 
resulting effective mass and the vector field strength are shown in Fig. |2|. One can see 
clearly that the effective mass decreases and the vector field strength increases rapidly as 



K increases. As shown in Ref. [^3|, the equivalent scalar and vector potentials appearing in 
the wave equation for a point-like nucleon are — M N and U v = g u uJ, respectively. Thus, 
our results indicate that the scalar and vector potentials for the nucleon become stronger as 
k gets larger. 

The corresponding in-medium bag radius is shown in Fig. [3[ In the QMC model, the 
bag radius decreases slightly in the medium. When k > 0, the bag constant drops relative 
to its free space value, which implies a decreasing bag pressure and hence gives rise to a 
increasing bag radius in the medium. With larger k, the bag radius grows more quickly. 
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FIG. 2. Results for the ratios M^/M^ and U v /Mn = g^TD/M^ as functions of the medium 
density, with Rq = 0.8 fm. The five curves correspond to, k = (solid), 7 (long-dashed), 8 
(dot-dashed), 9 (short-dashed), and 10 (dotted), respectively. The result from QHD-I is given by 
the open circles. 

This is consistent with the "swollen" nucleon picture drawn from the decrease of the meson 
masses in nuclear medium [|I7| -|26|| (see, however, Ref. ||27|| ). The in- medium bag constant 



is plotted in Fig. [|. One can see that when k > the bag constant decreases as density 
increases. The rate of this decrease gets larger for larger k values. Finally, the sensitivity 
of our results to the free space bag radius is illustrated in Fig. ||. For a given k value, one 
finds that the ratios B/Bq and M^/Mn increase and the ratio R/Ro and the vector field 
strength decrease as Ro increases. 

We observe that n ~ and k ~ 7 — 10 lead to qualitatively different physics. Unless 
one expresses the bag constant in terms of QCD operators and solves QCD in the nuclear 
matter, k is unknown. Nevertheless, one may get an estimate from the Brown-Rho scal- 
ing ansatz, $ ~ m* p /m p . Using m*/m p ~ 0.8 at = p% as suggested in Ref. and 



9 



1.8 
1.6 
1.4 

o 

K 

\ 
K 

1.2 
1.0 
0.8 

0.0 0.5 1.0 1.5 

FIG. 3. Result for the ratio R/Ro as a function of the medium density, with Rq = 0.8 fm. The 
five curves correspond to k = (solid), 7 (long-dashed), 8 (dot-dashed), 9 (short-dashed), and 10 
(dotted) , respectively. 

B/Bq ~ $ 4 , we get B/B ~ 0.4 at the saturation density. This requires n ~ 8.6, 9.1, 9.3 
for R = 0.6, 0.8, 1.0 fm, respectively. For these k values, the reduction of the bag constant 
dominates the attraction, the internal structure of the nucleon only plays a relatively minor 
role, and the large and canceling scalar and vector potentials for the nucleon appears nat- 
urally. Such potentials are comparable to those suggested by Dirac phenomenology [|],@], 
Brueckner calculations ||, and finite-density QCD sum rules Q, but somewhat smaller 
than those obtained in QHD-I. These large potentials also imply a strong nucleon spin-orbit 
potential. Therefore, we conclude that the essential features of relativistic nuclear phe- 
nomenology can be recovered when the decrease of the bag constant with increasing density 
is considered. 

The QMC model is valid only if the nucleon bags do not overlap significantly. With the 
k values suggested above, we find R/Ro ~ 1.25 at = p%- For R = 0.6 — 0.7 fm, as 
suggested by Guichon M, one finds R ~ 0.75 — 0.875 fm, which gives 47iR 3 p° N /3 ~ 0.3 — 0.48. 
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This indicates that the overlap between the bags is reasonably small at the saturation 
density, though a factor of two larger than in the usual QMC model. For larger R$ and/or 
higher densities, the overlap becomes more significant and the non-overlapping picture of 
the nuclear matter may become inadequate. However, it is unclear at this stage whether the 
overlap between the bags is effectively included in the scalar and vector mean fields. Further 
study is needed to clarify this issue. We also note that for the k values suggested above, 
the resulting nuclear incompressibility is comparable to that obtained in QHD-I, which is 
too large compared with the empirical value and that obtained in the usual QMC model. 
This may be fixed by introducing self-interactions of the scalar field, which, however, will 
introduce more free parameters. 

In summary, we have included the decrease of the bag constant in the quark-meson 
coupling model for the nuclear matter. This effectively introduces a new source of attraction, 
which needs to be compensated with additional vector field strength. When the change of 
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FIG. 4. Result for the ratio B/Bq as a function of the medium density, with Rq = 0.8 fm. The 
five curves correspond to k = (solid), 7 (long-dashed), 8 (dot-dashed), 9 (short-dashed), and 10 
(dotted) , respectively. 
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FIG. 5. Results for the ratios M^/M^ and U v /M^ = g^UJ/M^ as functions of the medium 
density, with k = 9.0. The three curves are for Rq = 0.6 fm (solid), 0.8 fm (dashed), and 1.0 fm 
(dot-dashed), respectively. 

the bag constant is large, as supported by partial chiral-symmetry restoration, large and 
canceling scalar and vector potentials for the nucleon emerge. The essential physics of the 
relativistic nuclear phenomenology can be recovered. The internal quark structure of the 
nucleon seems to play only a relatively minor role. On the other hand, the in-medium 
modification of the bag constant may play an important role in low- and medium-energy 
nuclear physics. The model presented in the present paper can be applied to variety of 
nuclear physics problems. Work along this direction is in progress and will be reported 
elsewhere. 
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